Distribution of voltage fluctuations in a current-biased conductor 
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We calculate the fluctuating voltage V(t) over a conductor driven out of equilibrium by a current 
source. This is the dual of the shot noise problem of current fluctuations I(t) in a voltage-biased 
circuit. In the single-channel case the distribution of the accumulated phase <E> = (e/h) JVdt is 
the Pascal (or binomial waiting-time) distribution — distinct from the binomial distribution of 
transferred charge Q — J Idt. The weak-coupling limit of a Poissonian -P(<£>) is reached in the limit 
of a ballistic conductor, while in the tunneling limit P($) has the chi-square form. 
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The current-voltage or charge-phase duality plays a 
central role in the theory of single-electron tunneling 
through tunnel junctions of small capacitance [[[]. At the 
two extremes one has a voltage-biased junction (in which 
the voltage is kept fixed by a source with zero internal 
resistance, while the current fluctuates) and a current- 
biased junction (fixed current from a source with infinite 
internal resistance, fluctuating voltage). In the current- 
biased case the Coulomb blockade induces periodic oscil- 
lations in the voltage Q , while in the voltage-biased case 
the Coulomb blockade is inoperative. 

Quantum mechanically, the duality appears because 
current I and voltage V are noncommuting operators || . 
This is conveniently expressed by the canonical commu- 
tator [$, Q] — ie of the transferred charge Q = f Q r I{t)dt 
and accumulated phase <i> = (e/h) V(t)dt (in a given 
detection time r). Moments of charge and phase deter- 
mine the measured correlators of current and voltage, 
respectively Q. 

While all moments of Q in a voltage-biased conductor 
are known (Q), the dual problem (moments of $ under 
current bias) has only been studied for the first two mo- 
ments (^,0. In the absence of Coulomb-blockade effects, 
the first two moments in the dual problems are simply 
related by rescaling I(t) —> V(t) x G (with G the con- 
ductance). One might surmise that this linear rescaling 
carries over to higher moments, so that the dual prob- 
lems are trivially related in the absence of the Coulomb 
blockade. However, the rescaling (as derived for example 
in Ref . || ) follows from a Langevin approach that is sus- 
pect for moments higher than the second — so that one 
might expect a more complex duality relation. 

The resolution of this issue is particularly urgent in 
view of recent proposals to measure the third cumulant 
of shot noise in a mesoscopic conductor §, [h], ||. Does 
it matter if the circuit is voltage biased or current biased, 
or can one relate one circuit to the other by a linear 
rescaling? That is the question addressed in this paper. 

We will demonstrate that, quite generally, the rescal- 
ing breaks down beyond the second moment. We cal- 
culate all moments of the phase (hence all correlators 
of the voltage) for the simplest case of a single-channel 



|l2| conductor (transmission probability T) in the (zero- 
temperature) shot-noise limit. In this case the charge 
Q = qe for voltage bias Vq = tufio/eT is known to have 
the binomial distribution Q 

p*M= Wpti-r)* 1 -' (i) 

We find that the dual distribution of phase $ = 2ircj) for 
current bias Jo = eqo /r is the Pascal distribution jl3| 

p <to (0 = (^ 1 1 )r»(i-r)*-*\ (2) 

(Both q and <fi & r e integers for integer 4>o and go-) 

In the more general case we have found that the distri- 
butions of charge and phase are related in a remarkably 
simple fashion for q,(f> — > oo: 

lnP,(0)=lnP («)+O(l). (3) 

(The remainder 0(1) equals \n(q/<f>) in the shot-noise 
limit.) This manifestation of charge-phase duality, valid 
with logarithmic accuracy, holds for any number of chan- 
nels and any model of the conductor. 

We have obtained these results in the course of an ex- 
tensive study of the quantum statistics of fluctuating cur- 
rents and voltages in electrical circuits containing meso- 
scopic conductors. Our method is the non-equilibrium 
Keldysh action technique. Within this technique, the 
electronic degrees of freedom are traced out and the 
quantum dynamics is represented by path integrals over 
fluctuating phases on the Keldysh contour Before 
presenting the derivation we give an intuitive physical 
interpretation. 

The binomial distribution (|l|) for voltage bias has the 
interpretation || that electrons hit the barrier with fre- 
quency eVo/h and are transmitted independently with 
probability T. For current bias the transmission rate is 
fixed at Io/e. Deviations due to the probabilistic nature 
of the transmission process are compensated for by an 
adjustment of the voltage drop over the barrier. If the 
transmission rate is too low, the voltage V(t) rises so 
that electrons hit the barrier with higher frequency. The 



number of transmission attempts ("trials") in a time r 
is given by (e/h) f£ V(t)dt = 4>. The statistics of the 
accumulated phase <f> is therefore given by the statistics 
of the number of trials needed for Iqt/c successful trans- 
mission events. This stochastic process has the Pascal 
distribution (g). 

Starting point of our derivation is a generalization to 
time-dependent bias voltage V(t) — (Ti/e)<b(t) of an ex- 
pression in the literature || [t4| for the generating func- 
tional Z[$(i), x{t)} °f current fluctuations: 



Z[9>,x\ - ^Texp|-y dt[$(t) + |x(*)]J(t) 
xTcxpj- f dt[-$(t) + §*(<)] 1(f) 



(4) 



(The notation T(T) denotes time-ordering of the ex- 
ponentials in ascending (descending) order.) Functional 
derivatives of the Keldysh action InZ with respect to 
x(t)/e produce irreducible correlators of the current op- 
erator I{t) to any order desired. This expression provides 
an example of the "tracing out" of degrees of freedom; 
The average is over electronic degrees of freedom, the 
electrical current being the only operator coupled to the 
phase. 

To make the transition from voltage to current bias 
we introduce a second conductor B in series with the 
mesoscopic conductor A (see Fig. Q). The generating 
functional Za+ b of current fluctuations in the circuit is 
a (path integral) convolution of Za and Zb, 

Z a +b{$,x} = J V^xVxiZa^x^Zb^-^X-Xi]- 

One can understand this expression as the average over 
fluctuating phases <&i, xi a t the node of the circuit shared 
by both conductors. 

In general the functional dependence of Za , Zb is 
rather complicated and non-local in time, but we have 
found an interesting and tractable low-frequency regime: 
The non-locality may be disregarded for sufficiently slow 
realizations of the fluctuating phases. In this regime the 
functional Z can be expressed in terms of a function S, 



lnZmt), X (t)} = / dtS(<Z>(t), X (t)) 



(6) 



The path integral (^|) can be taken in saddle-point ap- 
proximation, with the result 

S a+ b& x) = s A (i> s ,xs) + s B (i> - $ s , x - Xs). (7) 

Here $ s and Xs stand for the (generally complex) values 
of $i and xi at the saddle point (where the derivatives 
with respect to these phases vanish). 

Eqs. (6|) and (Q) are quite general and now we apply 
them to the specific circuit of Fig. [TJ. We assume that 
the mesoscopic conductor A (conductance G) is in sc- 
ries with a macroscopic conductor B with frequency de- 
pendent impedance Z(uS). (Any electromagnetic environ- 
ment with a Gaussian Keldysh action can be represented 
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FIG. 1: Mesoscopic conductor (shaded) in a circuit contain- 
ing a voltage source Vo and series impedance Z(ui). Both the 
current I through the circuit and the voltage drop V over 
the conductor may fluctuate in time. The dual problems con- 
trasted here are: voltage bias (Z — > 0, fixed V — Vo, fluctuat- 
ing 7) and current bias (Z — > oo, fixed / = Va/Z, fluctuating 
V). The phases \ appearing in Eq. (jB|) are indicated. 



by an external impedance.) The circuit is driven by a 
voltage source with voltage Vq. To avoid Coulomb block- 
ade effects Q, we assume that |Z(^)| <C h/e 2 at frequen- 
cies TkJ > max(e V, kT) (with T the temperature) . The 
zero-frequency impedance Z(0) = Zq = z$h/e 2 can have 
any value. 

Both the voltage drop V at the mesoscopic conductor 
and the current I through the conductor fluctuate in time 
for finite Z , with averages I = VqG(1 + ZqG) -1 , V — 
Vb(l + ZqG)^ 1 . Voltage bias corresponds to ZqG <C 1 
and current bias to Z G 3> 1, with I — Vq/Zq the im- 
posed current. There are three characteristic time scales: 
Tij max(eV r , kT), e/I, and the i?C-time of the circuit. 
(The inverse of the i?C-time is the frequency below which 
Z « Zq.) The low- frequency regime on which we concen- 
trate is reached for current and voltage fluctuations that 
are slow on any of these time scales. 

We assume that the temperature is sufficiently low 
(kT <C eV) to neglect thermal noise relative to shot 
noise. (Generalization to finite temperature is straight- 
forward, but will not be done here for simplicity.) The 
low-temperature, low- frequency Keldysh action of the ex- 
ternal impedance is simply S B (&,x) — ix^/^ z 0i while 
the action Sa of the mesoscopic conductor is given by || 



SA^,x) = ^-S(i X ), S(0 = £ln[l+(e*-l)T n ]. (8) 



n 
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The TVs are the transmission eigenvalues, with ~^2 n T n — 
Gh/e 2 = g the dimensionlcss conductance. (The formu- 
las are written without including spin or other degenera- 
cies @.) 

We seek the cumulant generating function of charge 



/ oo \ oo 

\g=0 / p=l P ' 



(9) 
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where ((q p )) is the p-th cumulant of the charge transferred 
during the time interval r. It is related to the Keldysh 
action (fjj) by 



F{t)=TS A+B (eV /h,-iO. 



(10) 



We also require the cumulant generating function of 
phase, Since V = Vo — ZqI (in the absence of 

thermal noise from the external impedance) , it is related 
to by a change of variables (from q to <p — (f>o —qzo). 
The relation is 



P =i p - 



(ii) 



In the limit Zq — > of voltage bias the saddle point 
of the Keldysh action is at $ s = <I>, \s = Xi an d 
from Ens. (0), (§), and (|ll|) one recovers the results 
of Ref. p|: The cumulant generating function J-oCO = 
tSa^Vq/H, — i£) = 0o<S>(£) and the corresponding prob- 
ability distribution 



1 d q 

p (q) = lim - — 
x^o q\ dx q 



N 



[J[l + (a:-l)r w ]' 



(12) 



The parameter 4>q = eVor/h is the number of attempted 
transmissions per channel, assumed to be an integer ^> 
1. The first few cumulants are (q)o = fog, ((q 2 ))o — 
<f>oJ2 n T n (l-T n ), ((q 3 )} = <£ £ n T„(l-T„)(l-2T„). In 
the single-channel case (N = 1, T\ = T) the distribution 
( |l2| ) has the binomial form (jl]). 

After these preparations we are now ready to generalize 
all of this to finite Zq, and in particular to derive the 
dual distribution of phase (0) under current bias. The 
key equation that allows us to do that follows directly 
from Eqs. (0) and @: 



(13) 



The implicit function cr(£) (which determines the saddle 
point of the Keldysh action) provides the cumulant gen- 
crating function of charge T for arbitrary series resistance 
z = (e 2 /h)Z . One readily checks that ,F(£) — > 4>oS(Q 
in the limit zq — > 0, as it should. 

By expanding Eq. ( |l3] ) in powers of £ we obtain a rela- 
tion between the cumulants {(q p )) of charge at Zq ^ and 
the cumulants {{q p ))o at Zq = 0. For example, to linear 
order we find (q) = (l+zog)~ 1 (q)o, which amounts to the 
obvious statement that the mean current / is rescaled by 
a factor 1 + ZqG as a result of the series resistance. The 
Langevin approach discussed in the introduction predicts 
that the same rescaling applies to the fluctuations. In- 
deed, to second order we find {(q 2 )) = (1 + zog)~ 3 ((q 2 ))o, 
in agreement with Ref. Q. 

However, if we go to higher cumulants we find that 
other terms appear, which can not be incorporated by 
any rescaling. For example, Eq. ( p"3| ) gives for the third 



cumulant 

«<z 3 » 



3z g (((g 2 ))o 

(1 + z g) 4 (1 + z g) 5 {q) 



«9 3 »c 



(14) 



The first term on the the right-hand-side has the ex- 
pected scaling form, but the second term does not. This 
is generic for p > 3: {(q p )) = (1 + zog)~' p ~ 1 ((q p )) plu s 
a non-linear (rational) function of lower cumulants p5[ . 
All terms are of the same order of magnitude in Zog, so 
one can not neglect the non-linear terms. 

Turning now to the limit zog — > oo of current bias, we 
see from Eq. (fL3j) that T — > with 



^oo(£)=<Zo£-<7oS inv (£Ao) 



(15) 



defined in terms of the functional inverse S lnv of S. The 
parameter qo — (f>a / zq = Ior/e (assumed to be an integer 
3> 1) is the number of charges transferred by the imposed 
current Iq in the detection time t. Transforming from 
charge to phase variables by means of Eq. ( |TT| ) , we find 
that Q — ► C?oo with 



<?oo(£) = -q S™(-0. 



(16) 



In the single-channel case Eq. ( |16| ) reduces to Goa(C) = 
—qo ln[l + r -1 (e~£ — 1)], corresponding to the Pascal dis- 
tribution (^). The first three cumulants are ((f)) = qo/T, 

U 2 )) = ( 90 /r 2 )(i - r), u 3 )) = (?o/r 3 )(i - r)(2 - r). 

For the general multi-channel case a simple expres- 
sion for P qo (cf>) can be obtained in the ballistic limit 
(all T n 's close to 1) and in the tunneling limit (all T„'s 
close to 0). In the ballistic limit one has Goo(0 = 
qo£,/N + qo(N — g)(e^' N — 1), corresponding to a Poisson 
distribution in the discrete variable N(f> — qo = 0, 1, 2, . . .. 
In the tunneling limit Goo(0 = — qo ln(l — £/<?), corre- 
sponding to a chi-square distribution P qo ((f)) oc <?o_1 e -9 * 
in the continuous variable <fi > 0. In contrast, the charge 
distribution P^ (g) is Poissonian both in the tunneling 
limit (in the variable q) and in the ballistic limit (in the 
variable N(f>$ — q) . 

For large q and (f), when the discreteness of these 
variables can be ignored, we may calculate P qo (4>) from 
Goo{£,) in saddle-point approximation. If we also calcu- 
late P(j> (q) from Fq(£) in the same approximation (valid 
for large (f>o and g), we find that the two distributions 
have a remarkably similar form: 

P^iq) = Ar 0o ( g )exp[rS(27r^o/T,g/r)] ! (17) 
P qo (<t>) = JV W) W)exp[rE(2^/r,<fo/T)]. (18) 



The same exponential function 



E(x, y) = Sa(x, -i£ s ) - y£s 



(19) 



appears in both distributions (with £ s the location of the 
saddle point). The pre-exponential functions and 
N qo are different, determined by the Gaussian integration 
around the saddle point. Since these two functions vary 
only algebraically, rather than exponentially, we conclude 
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FIG. 2: Comparison of the distributions of charge (dashed 
curve, with x — q/{q}) and of phase (solid curve, with x = 
tj>/{(j>)), calculated from Eqs. ( po[ ) and (^1) for AT = go = 
4>qT — 30 transferred charges in the tunne ing limit r < 1. 
The main plot emphasizes the non-Gaussian tails on a semi- 
logarithmic scale, the inset shows on a linear scale that the 
Gaussian body of the distributions coincides. 



that Eq. (||) holds with the remainder 0(1) = ]n{q/<j>) ob- 
tained by evaluating \n[2n{d 2 Y,/dx 2 ) 1 / 2 {d 2 Y,/dy 2 )- 1 /' 2 ] 
at x = 2ir(j)/T, y — q/r. 

The distributions of charge and phase are compared 
graphically in Fig. |2[ in the tunneling limit f < 1. We 



use the rescaled variable x = q/(q) for the charge and 
x = 4>/{4>) for the phase, and take the same mean number 
J\f = qa = 4>oT of transferred charges in both cases. We 
plot the asymptotic large- Af form of the distributions, 

Pchar g c(z) = (A^) 1 / 2 ^ 1 /^- 1 -* 111 ^, (20) 

PphascM = (AW^-leA^l-z+ln^ ( 2 1) 

corresponding to the Poisson and chi-square distribution, 
respectively. Since the first two moments are the same, 
the difference appears in the non-Gaussian tails. The 
difference should be readily visible as a factor of two in 
a measurement of the third cumulant: ((x 3 )) = Af~ 2 for 
the charge and ((x 3 )) = 2Af~ 2 for the phase. 

In summary, we have demonstrated theoretically that 
electrical noise becomes intrinsically different when the 
conductor is current biased rather than voltage biased. 
While the second cumulants can be related by a rescal- 
ing with the conductance, the third and higher cumulants 
can not. Experiments are typically carried out in an in- 
termediate regime, and for that purpose we have derived 
explicit expressions [such as Eq. (|l4[) ] for the crossover 
from voltage to current bias. The non-linear term that 
breaks the rescaling is of the same order of magnitude 
as the linear term, so it can not be neglected. From a 
fundamental point of view, the limit of full current bias is 
of particular interest. The counterpart of the celebrated 
binomial distribution of transferred charge Q turns out 
to be the Pascal distribution of phase increments. 

This work was supported by the Dutch Science Foun- 
dation NWO/FOM. 
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We record the result for the fourth cumulant, obtained by 
expansion of Eq. ( p^[ ) to order £ 4 : ({q 4 }} = (l + zo<?)~ 5 /M — 
10zo3(l + z g)' b fi2^a/fti + 15(z g) 2 (l + z g)~ 7 nl/n\, 
where we have abbreviated ({q p ))o = 



